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, Abstract We present the first Dyson- Sehwinger equation calculation of the light hadron spectrum 

that simultaneously correlates the masses of meson and baryon ground- and excited-states within a 
single framework. At the core of our analysis is a symmetry-preserving treatment of a vector-vector 
contact interaction. In comparison with relevant quantities the root-mean-square-relative-error/degree- 
of freedom is 13%. Notable amongst our results is agreement between the computed baryon masses and 
, the bare masses employed in modern dynamical coupled-channels models of pion-nucleon reactions. 

Our analysis provides insight into numerous aspects of baryon structure; e.g., relationships between 
the nuclcon and A masses and those of the dressed-quark and diquark correlations they contain. 
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1 Introduction 



(N 
> 

^SJ ■ Spectroscopy has long been a powerful means by which to expose the nature of. and interactions 

' between, the constituents of a compound system. In this context hadron spectroscopy has already 

produced many surprises. Most notable, perhaps, being the prediction of hadron substructure and its 
expression in terms of constituent-quark degrees-of-freedom, which were critical steps in the develop- 
ment of quantum chromodynamics (QCD). 

With QCD Nature has given us the sole known example of a strongly-interacting quantum field 
theory that is defined through degrees-of-freedom which cannot directly be detected. This empirical 
fact of confinement ensures that QCD is the most interesting and challenging piece of the Standard 
' Model. It means that building a bridge between QCD and the observed properties of hadrons is one 

^ , of the key problems for modern science. 
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In hadron physics the constituent-quark model has hitherto been the most widely applied spectro- 
scopic tool [J; and whilst its weaknesses are emphasized by critics and acknowledged by proponents, 
it is of continuing value because there is nothing better that is yet providing a bigger picture. On the 
other hand, steps are being taken with approaches that can rigorously be connected with QCD. For 
example, in the continuation of a more-than thirty- year effort, extensive resources are being invested 
in the application of numerical simulations of lattice-regularised QCD, with the claim that a spectrum 
which can reasonably be compared with experiment may soon be within reach 0; Q . Herein we bring 
a different, continuum perspective to computing the hadron spectrum; i.e., QCD's Dyson-Schwinger 
equations (DSEs) 

The DSEs have been applied extensively to the spectrum and interactions of mesons with masses 
less than 1 GeV [1; S • On this domain the rainbow-ladder approximation, which is the leading- 
order in a systematic and symmetry-preserving truncation scheme [1; Q, is an accurate and well- 
understood tool [m [S T13I: [13 that can readily be extended to explain properties of the light 
neutral pseudoscalar mesons [1^. Whilst the rainbow- ladder truncation is also reliable for ground- 
state heavy- heavy mesons (T6| . in order to make progress with heavy- light mesons it is necessary to 
employ the essentially nonperturbative truncation scheme introduced in Ref. (T7| . 

A Poincare covariant Faddeev equation was formulated in Ref. (l8l |. with a first exploration of 
its solution in Ref. [l9| . The equation is derived following upon the observation that an interaction 
which describes mesons also generates diquark correlations in the colour-3 channel [g^]. Despite the 
existence of this tractable truncation of the three-body problem in quantum field theory, the spectrum 
of baryons has almost escaped analysis via the DSEs. There are some notable exceptions; e. g., s tudies of 
the spectrum of ground-state octet and decuplet baryons using simple interaction kernels j2lt [22I : [23j , 
and an extension to baryons containing a single heavy-quark |24{ . An attempt has recently been 
made to unify the treatment of mesons and baryons through the consistent use of a rainbow-ladder 
truncation [25|. However, reliable studies are currently available only for the nucleon's ground-state. 
Our analysis is a modest step toward unifying a larger subset of meson and baryon spectra via a single 
interaction: whilst certainly not the last word, it should serve to illustrate the potential of the DSEs 
in this connection and provide reliable insights. 

In Sect[2]we formulate the bound-state problem for mesons and diquarks. Following Refs. [H;|13: we 
employ a vector- vector contact-interaction, which is regularised such that confinement is manifest, and 
chiral symmetry and the pattern by which it is broken are voraciously represented. Section [3] explains 
how the interaction's three parameters are constrained, and reports a computation of the masses for 
eight mesons and eight diquark correlations. Owing to a simplification we employ in constructing the 
Faddeev equation kernels, two additional parameters appear in our formulation of the nucleon and A 
Faddeev equations. This is described in Sect. SI which also details results for the masses of these states 
and their parity partners, and the radial excitations of those states. Section [S] provides a summary and 
perspective. 

2 Mesons and diquarks: formulating the bound-state problem 

The bound-state problem for hadrons characterised by two valence-fermions may be studied using the 
homogeneous Bcthe-Salpeter equation (BSE): 

[r{k; P)],„ = J ^[X{q; P)]srKi:{q, fc; P) , (1) 

where: F is the bound-state's Bethe-Salpeter amplitude and x(g; P) = S[q + P)rS{q) is its Bethe- 
Salpeter wave-function; r,s,t,u represent colour, flavour and spinor indices; and K is the relevant 
fermion-fermion scattering kernel. This equation possesses solutions on that discrete set of P^-values 
for which bound-states exist. (Our Euclidean metric conventions are specified in App.Kl) In Eq. 
S is the dressed-quark propagator; viz.. the solution of the gap equation: 

j^^9^D^^{p- q)—-i^S{q) — r^{q,p), (2) 

wherein m is the Lagrangian current-quark mass, D^j/ is the gluon propagator and is the quark-gluon 
vertex. 
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2.1 Rainbow-ladder truncation 

For ground-state, charged pseudoscalar- and vector-mesons constituted from a valence-quark and 
-antiquark with equal current-mass, the rainbow-ladder truncation of the Bethe-Salpeter and gap 
equations provides a good approximation [HI; This means r^{p,q) = in both Eq. ^ and 
the construction of K in Eq. ([T]) , so that one works with 

/d^q A" A" 

■J^^9^D^,^{p-q) — -i^S{q) — -iy{q,p), (3) 

nk;P) = - J^g'D,Ap-q)Y^,S{q + P)riq;P)Siq)^j,. (4) 

In this truncation, colour- antitriplet quark-quark correlations (diquarks) are described by an ho- 
mogeneous Bethe-Salpeter equation that is readily inferred from Eq. viz. (20| . 



r,,{k-P)H'^ = - j ^g^D^,ip-q)^j^S{q + P)r,,iq;P)H^[S{-q)f 



4^ \a \' \a~\^ 

b.f , (5) 



i2n) 

where c = 1, 2, 3 is a colour label and {-ff^} are defined in Eq. (jC.6[) . Using the properties of the Dirac 
and GcU-Mann matrices, it is straightforward to show that 

r,,ik;P)C^ = -\J ^g^D,,Ap-q)^j,,S{q + P)r,,iq;P)C^Siq)^J,, (6) 

which explicates the observation made in the Introduction; i.e., an interaction that binds mesons also 
generates strong diquark correlations in the colour-3 channel. It follows moreover that one may obtain 
the mass and Bethe-Salpeter amplitude for a diquark with spin-parity from the equation for a 
J^^-meson in which the only change is a halving of the interaction strength. The flipping of the sign 
in parity occurs because fermions and antifermions have opposite parity. 

We note that the rainbow-ladder truncation usually generates asymptotic diquark states. Such 
states are not observed and their appearance is an artefact of the truncation. Higher-order terms in 
the quark-quark scattering kernel, whose analogue in the quark-antiquark channel do not materially 
affect the properties of vector and flavour non-singlet pseudoscalar mesons, ensure that QCD's quark- 
uark scattering matrix does not exhibit singularities which correspond to asymptotic diquark states 
IE Ell- Nevertheless, studies with kernels that don't produce diquark bound states, do support a 
physical interpretation of the masses, "T-(gg) p , obtained using the rainbow-ladder truncation; viz., the 
quantity i{qq) p ■= l/m(gg) ^ may be interpreted as a range over which the diquark correlation can 
propagate before fragmentation. 



2.2 Vector-vector contact interaction 

References [26t [27| have shown that a momentum-independent interaction of vector x vector character 
is capable of providing a description of tt- and p-meson static properties which is comparable to that 
obtained using more-sophisticated, QCD-renormalisation-group-improved interactions [l3;[2^;[2^. This 
is sufficient justification for proceeding with an analysis of mesons and diquarks using 

g'^Df,^{p -q)^ (5^1^^ , (7) 

where mc is a gluon mass-scale. (Such a scale is generated dynamically in QCD [sO; [13; [H; [s^-) 
With this interaction the gap equation becomes 

^''^^='^-P + '^+t^J^^^^ 5(g) 7, . (8) 
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The integral possesses a quadratic divergence, even in the chiral Umit. If the divergence is regularised 
in a Poincare covariant manner, then the solution is 



(9) 



where M is momentum-independent and determined by 

M , 1 



M = 771 



s + NP 



(10) 



To continue, one must specify a regularisation procedure. We write [3£ 



M2 



M2 



where Tji-^uv are, respectively, infrared and ultraviolet regulators. It is apparent from Eq. that a 
nonzero value of Tir =: 1/vlir implements confinement by ensuring the absence of quark production 
thresholds [H;|33|- Furthermore, since Eq. (O does not define a rcnormalisable theory, A^-^, := I/tuv 
cannot be removed but instead plays a dynamical role and sets the scale of all dimensioned quantities. 
The gap equation can now be written {r(a,y) is the incomplete gamma-function) 



M = 771 



M 



C'"(A/') , C'"(Af^) = A/'[r(-l, Af-^T^J - r(-l, A/Vj^.)] 



(11) 



Using the interaction we've specified, the homogeneous BSE for a pseudoscalar meson is 

4 1 f d^q 



S 771 



G 



(12) 



With a symmetry-preserving regularisation of the interaction in Eq. ([7]), the Bethe-Salpeter amplitude 
cannot depend on relative momentum and hence may be written 



^o-(^)=75 



(13) 



Crucially, the amplitude contains Fq-{P), a part of pseudovector origin. It is an essential component 
of a pseudoscalar meson, which has significant measurable consequences [26t [27l : l38l | and thus cannot 
be neglected. 

Following the discussion in Sect. 12. 11 it is straightforward to write the Bethe-Salpeter equation for 
a = diquark; viz., 



r: 



99„+^ > 3„j2 



G 



(27r)4 



^^S{q + P)r^^^^{P)S{q)^^, 



(14) 



where 



(15) 
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2.3 Ward-Takahashi identities 



In studies of the hadron spectrum it is critical that a computational approach satisfy the vector and 
axial-vector Ward-Takahashi identities. Without this it is impossible to preserve the pattern of chiral 
symmetry breaking in QCD and hence a veracious understanding of hadron mass splittings is not 
achievable. The 771 = axial-vector identity states (fc+ = k + P) 



P,,r5^{k+, k) = 5-i(fc+)i75 + ii^s-\k) 
where /5^(fc+, k) is the axial- vector vertex, which is determined by 



(16) 



4 1 

Oil Lf^ 



——^jaS{q+)r5f,{q+,q)S{q)ja 



(17) 



One must therefore implement a regularisation of this inhomogeneous BSE that maintains Eq. (|16[). 

To see what this entails, contract Eq. ([T7)) with and use Eq. (jl6p within the integrand. This 
yields the following two chiral limit identities: 



M - - — 
3 



= 



dS 
(2^ 



(2^)4 



1 



AP q\ 
P-q 



AP 



M2 



(18) 
(19) 



which must be satisfied after regularisation. Analysing the integrands using a Feynman parametrisation, 
one arrives at the follow identities for P^ = = m: 



A4 




16 M 

3 TTlp 



d'^q 



1 



dS W 



(2^)4 [g2 + M2] ' 

AP 



(27r)4 [92 + APf 



(20) 
(21) 



Equation (j20p is just the chiral-limit gap equation. Hence it requires nothing new of the regular- 
isation scheme. On the other hand, Eq. (|2ip states that the axial-vector Ward-Takahashi identity is 
satisfied if, and only if, the model is regularised so as to ensure there are no quadratic or logarithmic 
divergences. Unsurprisingly, these are the just the circumstances under which a shift in integration 
variables is permitted, an operation required in order to prove Eq. (|16p . 

We observe in addition that Eq. (|16p is valid for arbitrary P. In fact its corollary, Eq. ([T5)) . can be 
used to demonstrate that in the chiral limit the two-flavour scalar-meson rainbow-ladder truncation 
of the contact-interaction DSEs produces a bound-state with mass TOo- ~2Ad [s^ (see Apd. IB.31) . The 
second corollary, Eq. (|19p , entails 



0^ da [C'"(u;(Af2^a,p2)) + C['\ujiM^ ,a, P^))] , 
Jo 



(22) 



with cj(A/^a,p2) ^ aP + a{l - a)P'^ and Ci"(z) = -z(d/dz)C™(z). 



3 Mesons and diquarks: computed masses 

In App.[B]we present explicit forms of the homogeneous rainbow-ladder Bethe-Salpeter equations for 
ground-state = 0~ , 0+, 1^, 1+ mesons and diquarks, obtained using the interaction and regularisa- 
tion scheme described above. In order to present numerical results, the values of our parameters must 
be fixed. 
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Table 1 Results obtained with (in GeV) ma ~ 0.132, Air = 0.24, yluv = 0.905, which yield a root-mean- 
square relative-error of 13% in comparison with our specified goals for the observables. Dimensioned quantities 
are listed in GeV. 



m 


En 


Fn 


Ep 


M 






mp 




fp 





3.568 


0.459 


1.520 


0.358 


0.241 





0.919 


0.100 


0.130 


0.007 


3.639 


0.481 


1.531 


0.368 


0.243 


0.140 


0.928 


0.101 


0.129 



Table 2 Meson masses (GeV) computed using our contact-interaction DSE kernel, which produces a 
momentum-independent dressed-quark mass M = 0.37 GeV from a current-quark mass of m = 7MeV. "RL" 
denotes rainbow-ladder truncation. Row-2 is obtained by augmenting the RL kernel with repulsion generated 
by vertex dressing (see Eq. (|25[) and associated text). The text ar ound Eq. (|30[) explains the errors on the 
masses of radially-excited states. Row-3 lists experimental masses [l^ for comparison. NB. We implement 



isospin symmetry 


so, e.g 




— TUp, m/j = 




etc. 














rrin 


nip 


771 o- 






m,r* 


TTlp. 




a, 


RL 

RL * gio 




0.14 
0.14 


0.93 
0.93 


0.74 
1.29 


1.08 
1.38 




1.38 ±0.06 
1.38 ±0.06 


1.29 ±0.07 
1.29 ±0.07 


1.41 ±0.06 
1.47 ±0.04 


1.30 ±0.06 
1.47 ±0.03 


experiment 




0.14 


0.78 


0.4 - 1.2 


1.24 




1.3±0.1 


1.47 


1.2 - 1.5 


1.43 



3.1 Ground states 

The interaction defined in Sect. 12. 21 possesses three parameters: Air, ^uv, and tug- We fix ylir = 
0.24 GeV w ylqcD, since rqcD := 1/^qcd ~ 0.8 fm is a length-scale typical of confinement, so that 
only the other two parameters are active. We fix their values by performing a least-squares fit in the 
chiral limit to M° = 0.40 GeV, m° = 0.78 GeV, f° = 0.088 GeV, /° = 0.15 GeV k° = (0.22 GeV)^. 
The hitherto undefined entries in this list are: the light-meson leptonic decay constantsQ 

f. = ^^[E^~2K]]C-MP'^~rnl), f^^J-^KP{~m^) (23) 

and the in-pion condensate [4l| 

= I.^[E.1CIe{-^iI) + IClpi^ml)] . (24) 

In these expressions, iJ^, Ft, and Ep are the canonically- normalised Bethe-Salpeter amplitudes (see 
App.[B]). This procedure yields the results in Table[Tl and the masses of the meson and diquark ground 
states reported in Tables [H and 121 

The pattern of meson masses is typical of the rainbow- ladder truncation (28l : [291 : |43| : tt- and p- 
mesons are described well but their parity partners - the a- and ai-mcsons - are not. The origin 
and solution of this longstanding puzzle are now available following a novel reformulation of the BSE 
[Trj . which is valid and tractable when the quark-gluon vertex is fully dressed. In employing this 
approach to study the meson spectrum it was found that DCSB generates a large dressed-quark 
anomalous chromomagnetic moment and consequently that spin-orbit splitting between ground-state 
mesons is dramatically enhanced [4^ : l45l | . This is the mechanism responsible for a magnified splitting 
between parity partners; namely, essentially-nonperturbative DCSB corrections to the rainbow-ladder 
truncation largely-cancel in the pseudoscalar and vector channels but add constructively in the scalar 
and axial- vector channels. 

With this in mind, we introduced spin-orbit repulsion into the scalar- and pseudovector-meson chan- 
nels through the artifice of a phenomenological coupling (7gQ < 1, introduced as a factor multiplying 
the kernels defined in Eqs. (|RT9)) . (|R26l) . The valucO 

gso = 0.240 (25) 

^ These expressions may be computed in a straightforward manner from the general formulae in, e.g., Ref. [43| . 
^ NB. gso ~ 1 means no repulsion. The mass changes slowly with diminishing gso\ e.g., gso ~ 0.50 yields 
mai= 1.23 GeV. 
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Table 3 Diquark masses (GeV) computed using our contact-interaction DSE kernel, which produces a 
momentum-independent dressed-quark mass A'l = 0.37 GeV from a current-quark mass of m = 7MeV. "RL" 
denotes rainbow-ladder truncation. Row-2 is obtained by augmenting the RL kernel with repulsion generated 
by vertex dressing (see Eq. H25|) and associated text). 









mqqa_ 








"^^^^ 




RL 

RL * glo 


0.78 
0.78 


1.06 
1.06 


0.93 
1.37 


1.16 
1.45 


1.39 ±0.06 
1.39 ±0.06 


1.32 ±0.05 
1.32 ±0.05 


1.42 ±0.05 
1.50 ±0.03 


1.33 ±0.05 
1.52 ±0.02 



is chosen so as to obtain the experimental value for the ai-p mass-splitting, which we know to be 
achieved by the corrections described above [ij 0; . This expedient produces the results in Row-2 
of Tables [5] and [3] It is noteworthy that the shift in ruai accompanied by an increase of TOo- and that 
the new value matches an estimate for the gg-componcnt of the cr-meson obtained using unitarised 
chiral perturbation theory [46l |. 

Tensor mesons are made conspicuous by their absence from Table [21 This is readily explained. In 
constituent-quark models one may only construct a J = 2 state constituted from two J = ^ quarks if 
the system's ground state contains at least one unit of orbital angular momentum. The analogue of this 
statement in quantum field theory is expressed in the requirement that a normal tensor meson's Bethe- 
Salpeter amplitude must depend at least linearly on the relative momentum [47^ . This is impossible 
using a symmetry-preserving regularisation of the interaction in Eq. ([7]) and hence it doesn't generate 
tensor meson bound states. 

Some remarks on the spectrum of ground-state diquarks are also appropriate here. Our computed 
values for the masses of the scalar and axial-vector diquarks are commensurate with other estimates 
based on the rainbow-ladder truncation [3|4^; and with numerical simulations of lattice-QCD js^]. In 
addition, it is noteworthy that our results for the gso-corrected masses of the ground-state pseudoscalar 
and vector diquarks are just 10% smaller than the values determined in Ref. [l^. Until recently the 
separable model employed therein for the Bethe-Salpeter kernel was unique in providing a realistic 
value for the ai-p mass-splitting [sij . 

It is also of interest to elucidate the role of the scalar diquark's vector component; i.e., the Fqq^^- 
term in its Bethe-Salpeter amplitude. Absent this term, the diquark mass drops by 62 MeV. Its presence 
therefore produces a small amount of repulsion. In Ref. [4^ this vector-component of the Bethe-Salpeter 
amplitude was found to provide a repulsive shift of 83 MeV in the scalar diquark's mass, whilst in 
Ref. m the shift is +80 MeV. 



3.2 Meson and diquark radial excitations 

In quantum mechanics the radial wave function for a bound-state's first radial excitation possesses a 
single zero. A similar feature is expressed in quantum field theory: namely, in a fully covariant approach 
a single zero is seen in the relative-momentum dependence of the leading Tchebychev moment of the 
dominant Dirac structure in the bound state amplitude for a meson's first radial excitation [52l |. 
The existence of radial excitations is therefore very obvious evidence against the possibility that the 
interaction between quarks is momentum- independent: a bound-state amplitude that is independent 
of the relative momentum cannot exhibit a single zero. One may also express this differently; namely, 
if the location of the zero is at /cq, then a momentum-independent interaction can only produce reliable 
results for phenomena that probe momentum scales <C kg. In QCD. feg ^ M and hence this criterion 
is equivalent to that noted in Ref. [1^. 

Herein, however, we skirt this difficulty by means of an expedient employed in Ref. [ssj : i.e., we 
insert a zero by hand into the kernels defined in Eqs. (jB.2|) - (jB.5() . (|B.14p . (|B.19() . (|B.26|) . This means 
that we identify the BSE for a radial excitation as the form of Eq. Q obtained with Eq. ([7]) and 
insertion into the integrand of a factor 

1 - d^q^ , (26) 

which forces a zero into the kernel at = where djr is a parameter. It is plain that the presence 

of this zero has the effect of reducing the coupling in the BSE and hence it increases the bound-state's 
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mass. Although this may not be as transparent with a more sophisticated interaction, a quahtatively 
equivalent mechanism is always responsible for the elevated values of the masses of radial excitations. 

To illustrate our procedure, consider the BSE for the vector meson, in which the following replace- 
ment is made 13 

Kf{P^)^K'''{P^)^-^ f daa{\~a)P^T^{uj{M\a,P^)), (27) 

6lT TJIq Jq 

where 

^'"(w(M2,a,F2)) =C'"(w(M2,a,p2)) _rf^pi"(^(Af2^a,p2)), (28) 
V^^{Lo{M-\a,P^))= dss^-^jj^^ j ^ dT-c^^[-TUj{M\a,P^)\, (29) 

r^{z) ^ ^z{d/dz)r''{z) and 7i(z) = Ti{z)/z. 

Regarding the location of the zero, motivated by extant studies in the pseudoscalar channel [52l |. 
we choose 1/djr = M^. The position of the zero in the leading Tchcbychev moment of an excited state 
in a given channel is an indication of that state's size. Hence it is an oversimplification to place the 
zero at the same location in each channel. Therefore, in Tables [2] and [Sj we report results obtained 
subject to a 20% variation in the zero's location; i.e., determined with 

-!- =Af2(-i o±0.2). (30) 
dj^ 

In order to evaluate the credibility of masses we compute subsequently for baryons, it is important 
to consider critically the results in the right columns of Table [2l In this connection it is noteworthy that 
whilst no parameters were tuned the computed masses are in good agreement with the known spectrum. 
Hence, this simple model produces a phenomenology which represents a considerable improvement over 
that of existing DSE studies. We therefore judge it to provide a solid foundation for a study of baryons. 



4 Spectrum of Baryons 

We compute the masses of light-quark baryons using a Faddeev equation built from the interaction in 
Eq. ([7]) and the diquark correlations discussed quantitatively in Sect. [31 The general structure of this 
equation for nucleon and A states is described in App. lC.ll 

The bound-state equations specific to our model are defined once the detailed forms of the kernels 
in Eqs. (|C.16|) and (jC.22l) arc specified; namely, by the structure of the dressed-quark propagator, 
the diquark Bethe-Salpeter amplitudes and the diquark propagators. In completing these kernels we 
make a drastic simplification; viz., in the Faddeev equation for a baryon of type B = N, A, the quark 
exchanged between the diquarks is represented as 

(31) 

where qb is discussed below. This is a variant of the so-called "static approximation," which itself 
was introduced in Ref. ^3] a-nd has subsequently been used in studies of a range of nucleon proper- 
ties (ssf . In combination with diquark correlations generated by Eq. ([7]), whose Bethe-Salpeter ampli- 
tudes are momentum- independent, Eq. pip generates Faddeev equation kernels which themselves are 
momentum-independent. The dramatic simplifications which this produces are the merit of Eq. (j31|) . 



^ The procedure is a little more involved in the pseudoscalar channel owing to the axial- ve ctor Ward- 
Takahashi identity and the richer structure of the Bethe-Salpeter amplitude. It is detailed in App. lB.ll 



9 



4.1 Ground-state A and nucleon 

Owing to its inherent simplicity, we use the A to ihustrate the construction of a Faddeev equation. 
With a momentum-independent kernel, the Faddeev amphtude cannot depend on relative momentum. 
Hence Eq. (|C.15|) becomes 

V,p{e; P)up{P) = f^{P) /d u,{P) . (32) 

NB. Regarding Eq. (|C.15() in general, one might naively suppose that isospin-one tensor diquarks could 
play a material role in the Faddeev amplitude for a ground state A. However, this notion can quickly be 
discarded because ground-states are distinguished by containing the smallest amount of quark orbital 
angular momentum, L, and a tensor diquark is characterised by L > 1. 
Using Eq. Eq. (|a22]) can be written 



r{P) u,{P) = 4^ y ^ M^,{t, P) r{P) u^P) , (33) 

with {K = -e + P, P^ = -m^) 

M^M;P) = 2^^l\K)^^l\-P)S{^)Af,{K) , (34) 

where the "2" has arisen through the isospin contractions. 

At this point, one post-multiplies by up{P]r) and sums over the polarisation index to obtain, 
Ea. (|XT3l) . 



yl+(F)i?^,(P)=4^ j ALm^,{1,P)A+{P)R,,{P)., 
which, after contracting with S^p, yields 



(35) 



1 = tro y ^ M^A^; P) A+{P)R,,iP) (36) 



8 9a f I f-f.Pl'i'm^ fn^+iK-P)"^] 

+mA[2m.At ■ KK ■ P + ZM{ml^^^m\ + {K ■ Pf)]^ , (37) 



where Eqq_^^{K) is the canonically-normalised axial-vector diquark Bethe-Salpeter amplitude, 
Eq. (|B.16p . Now, with the aid of a Feynman parametrisation, the right hand side becomes 



8 9l f dH 



d^ TU ^^7 ^ -^ P [3 m'm'^ + {K ■ P)'] 



3Mm^m2^^^y (271)4^0 [{t - aPy + aA{a,M,mqq^^,mA)Y \ ' ""^^ 

^mA[2mAi ■KK-P + MI{mlg^^ml^ + {K ■ Pf)]^ (38) 

where 

cr/i(a. A/, mqq^^ , itla) = (1 — a) il/^ + aniqq^^ — a{l — a) rn^. (39) 

We employ a symmetry-preserving regularisation scheme. Hence the shift £ — >■= £'+aP is permitted, 
whereafter 0(4)-invariance entails • P = so that one may set 

e-P^aP^ , A' • P = (1 - a)p2 ^ e-K ^a{l- a)P^ , (40) 

and therewith obtain [56j 

M m^q^^ J i2ny Jo " P + ^^(a, Af, m,,^^ , mzi)]^ ^ ' 

= 17^^Tt/ da{ml +{l-afml){amA+M)C"{aA{a,M,mqq^^,mA)). (42) 

m ITlqq^^ ZTT Jq 
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Fig. 1 Solid curve, 5m := {niA — raN); dashed 
curve, ruA', dot-dashed curve, mjv - all plotted as a 
function of Sqq^^g := (m^g^^ — mqq^^). The vertical 
dotted line marks the model's predicted values when 
qn = 1 = QA- 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 

mi+ - mo+(GeV) 

This is an eigenvalue problem whose solution yields the mass for the dressed-quark-core of the A- 
resonance. If one sets = 1, then rriA = 1.60 GeV. 

Construction of the explicit form for the nucleon's Faddeev equation is a straightforward general- 
isation of the procedure used above for the A. However, it is algebraically more complicated and in 
App. lC.2l we simply present the result. If one sets qm ~ 1, then m^r = 1.27 GeV. It is noteworthy that 
with Fqq^^ = 0, one obtains tojv = 1.14. This comparison shows that the scalar-diquark's vector- like 
component produces approximately 130 MeV of repulsion within the nucleon. 

One may now read that in the truncation we've described thus far 

5m := {rriA ~ niN) — 0.33 GeV cf. Sqq^g {mqq^^ — mqq^^) — 0.28 GeV. (43) 

We note that these mass differences are correlated; and both vanish together in the limit of infinitely 
heavy current-quark masses, approaching zero from above as the current-quark mass increases [STj . The 
latter results are a model-independent consequence of heavy-quark symmetry, kindred to the behaviour 
of the mass-splitting between vector and pseudoscalar mesons [l^. 

The causal connection between 5m and (^gg^j, is readily illustrated. For example, with all other 
elements held fixed, the latter determines the former. This is made plain in Fig.[Tl which depicts the 
mass-difference 5m = {mA — m^) and the masses mA, m^, all as a function of Sqq-^^, which was 
reduced from its model-preferred value by increasing the coupling in the BSE for the axial-vector di- 
quark whilst keeping all other couplings and masses constant. The behaviour of each curve is readily 
understood. The A is an uncomplicated bound state composed of a dressed-quark and an axial-vector 
diquark. Since diquark breakup and reformation mediated by the exchange of a dressed-quark is at- 
tractive, then decreasing the mass of the axial- vector diquark increases the amount of attraction in the 
channel because the attraction operates over longer range. There is nothing in this channel with which 
the increased attraction can compete, hence mA drops rapidly with decreasing 6qq-^g. The nucleon is 
more complicated. Its kernel expresses interference between quark-exchange in the scalar and axial- 
vector diquark channels, which provides resistance to change because the scalar-diquark properties are 
held fixed. Hence, m^ drops more slowly with decreasing Sqq-^g. Therefore the behaviour of the mass 
difference 5m is driven primarily by the change in mA- 



4.2 Pion loops, m^ and mA 

The results described hitherto suggest that whilst corrections to our truncated DSE kernels may have a 
material impact on m^r and mA separately, the modification of each is approximately the same, so that 
the mass-difference, 5m, is largely unaffected by such corrections. Indeed, this is consistent with an 
analysis |58| that considers the effect of pion loops, which are explicitly excluded in the rainbow-ladder 
truncation [1^: whilst the individual masses are reduced by roughly 300 MeV, the mass difference, 5m, 
increases by only 50 MeV. 

We emphasise that it is essential not to miscount when incorporating the effect of pseudoscalar 
meson loops. In practical calculations these effects divide into two distinct types. The first is within 
the gap equation, where pseudoscalar meson loop corrections to the dressed-quark-gluon vertex act 
to reduce uniformly the mass- function of a dressed-quark [l^ [s^ ■ This effect can be pictured as a 
single quark emitting and reabsorbing a pseudoscalar meson. It can be mocked-up by simply choosing 
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m^(GeV2) m^(GeV2) 

Fig. 2 Left panel - Evolution with current-quark mass of the: nucleon mass, mjv (solid curve); scalar-diquark- 

only nucleon mass, m'jl^ o'^ly (dashed curve); and A mass, tua- Right panel - Ratio ^^^■^/miv as a function 
of current-quark mass. It drops by just 5% on the domain depicted. In all panels here and below the current- 
quark mass is expressed through the computed value of m^; = 0.49 GeV^ marks the s-quark current- mass. 



the parameters in the gap equation's kernel so as to obtain a dressed-quark mass-funetion that is 
charaeterised by a mass-scale of approximately 400 MeV. Such an approach has implicitly been widely 
employed with phenomenological success [1; [a; j ^^'^ that which we employ herein. 

The second type of correction arises in connection with bound-states and may be likened to adding 
pseudoscalar meson exchange between dressed-quarks within the bound-state 5^; 60l : l6ll : [6^ . as opposed 



to the first type of effect; i.e., emission and absorption of a meson by the same quark. The type-2 
contribution is that computed in typical evaluations of meson-loop corrections to hadron observables 
based on a point-hadron Lagrangian. This fact should be borne in mind when using formulae, such 
as those in Ref. [SSf . to estimate the size of meson- loop corrections to bound-state masses computed 
using the DSEs. 

This discussion establishes that it is correct to use such formulae herein, just as it was in Refs. [6^ 
1131 • Their straightforward application using a common meson-baryon form- factor mass-scale of 0.8 GeV 
yields a shift of (—300 MeV) in m^v and (—270 MeV) in m/i, from which one may infer that our type-2 
uncorrected Faddeev equations should produce rriN = 1.24 GeV and m/i = 1.50 GeV, values which are 
plainly of the appropriate size. For the Z\-resonance there is another estimate, which is arguably more 
sophisticated. Namely, that produced by the Excited Baryon Analysis Center (EBAC) [65[, which used 
a realistic coupled-channels model to remove meson dressing from the A and expose a dressed-quark- 
core bare-mass of 1.39 GeV. Following these observations we return to Eq. (pij) and choose 

5Ar = 1.18, .9zi = 1-56 ^ = 1.14 GeV, mzi = 1.39 GeV, (5m = 0.25 GeV (44) 

because the listed outcomes of this choice are consistent with the information presented above and 
Refs. [6i|:|6l:|6l. 



4.3 Evolution of ground-state masses with current-quark mass 

Within a framework such as we employ it is straightforward to map the evolution with increasing 
current-quark mass, m, of bound-state masses and the splittings between them. In rainbow-ladder 
truncation it is only the current-mass which varies because the dressed-quark-gluon vertex is indepen- 
dent of m. Notably, in systems related to the pseudoscalar- and vector-meson channels, non-resonant 
corrections to the rainbow-ladder truncation do not materially affect mass-splittings [ll[. Hence the 
interpretation of our results on splittings between the bound-states' quark-cores should be robust. 

The nucleon is constituted from scalar and axial-vector diquark correlations. It is therefore of 
interest to determine the effect of the axial-vector correlation on the nucleon's mass. This can be 
read from Fig.[5J the axial- vector diquark-correlation provides attraction in the nucleon channel, in an 
amount which is almost independent of current-quark mass up to values matching the s-quark current- 
mass. This attraction ensures the nucleon remains lighter than the A for all values of current-quark 
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^(GeV^) mf(GeV2) 

Fig. 3 Left panel - Evolution with current-quark mass of the: nucleon mass, mjv {solid curve); the sum 
[M + rUgq^^] {dashed curve); and 3M {dotted curve). Right panel - Evolution with current-quark mass of the 
ratio mjv/[3Af], which varies by less-than 1% on the domain depicted. 




Fig. 4 Left panel - Evolution with current-quark mass of the: A mass, tua {solid curve); and [AL + nigq ^] 
(dashed curve). Right panel - Evolution with current-quark mass of the ratio mA/[M + mqq_^^], which does 
not vary noticeably on the domain depicted. 

mass. A nucleon constituted solely from a dressed-quark and scalar-diquark correlation will finally 
become more massive than the A. In our case, this occurs for m > 0.8ms. 

In the left panel of Fig. [3] we depict the evolution with current-quark mass of toat, the mass of the 
nucleon's dressed-quark-core, and compare it with the evolution of the combinations [Af -t- Triqq^^ ] and 
3M. In conjunction with the ratio in the right panel, it is evident that ttiat is given by 3Af plus a 
small contribution that is almost independent of current-quark mass0 Hence, if one inflates M in an 
attempt to anticipate type-1 pseudoscalar-meson vertex-corrections within the gap equation [l3j . then 
the nucleon mass will increase commensurately. It is noteworthy that for m > rUs, fnqq^j^ ~ 2M . 

In the left panel of Fig.|3]we display the evolution with current-quark mass of m^, the dressed- 
quark-core mass of the Z\, along with that of the sum M + 'm-qq-^+ ■ As the right panel makes plain, to 
a very good level of approximation = AI + mqq_^^ . Recall now that rnqq_^^ > rnqq^^ for all values 
of current-quark mass, with [mqq_^^ — "niqq^^] — >■ 0^ as m — >■ oo. It follows that if one inflates M in an 
attempt to anticipate type-1 pseudoscalar-meson vertex-corrections within the gap equation [l3|, then 
the A mass will increase by an amount far larger than is seen with the nucleon. Indeed, the amount will 
be commensurate with the inflation in 3A/ plus the scaled increase of the bound-state energy-excess 
— ] • This must be understood if erroneous conclusions about the nature of the quark-core 

of the A are to be avoided. 



In this near-proportionality there is a similarity with models of the constituent-quark type. NB. It is our 
confining regularisation of the contact-interaction which enables the mass of the bound-state to be greater than 
that contained in the masses of its constituents. 
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Fig. 5 Left panel ~ Evolution with current-quark mass of the: zi-nucleon mass difference [tua — mjv] (solid 
curve); and axial-vector-scalar-diquark mass difference [mqq_^_^ ~''^qqQ+] (dashed curve). Right panel - Evolution 
with current-quark mass of the ratio [m^ — m]v]/[mqq^^ — rriqq^^], which falls by 40% on the domain depicted. 



28 r 



Fig. 6 Evolution with current-quark mass of the 
contribution made by axial-vector-diquark correla- 
tions to the unit-normalisation of the nucleon's Fad- 
deev amplitude. At the model-preferred current- 
quark mass, the probability is 22%. This rises to 
26% at m = nia. 
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Figure [S] depicts the evolution with current-quark mass of the mass differences Sm = [rriA — mjv] 
and 5<j5ii = ~™ij«o+]- seen elsewhere [57| . with increasing m, Sm decreases more rapidly 

than Jqq^i and both mass-differences approach zero uniformly from above. The difference between 
the rates of decrease is readily explained. The difference 6qq^^ becomes smaller because increasing the 
quark mass suppresses hyperfine interactions, as demonstrated elsewhere [l^ for the analogous case 
of the vcctor-pseudoscalar-meson mass difference. This reduction introduces circumstances analogous 
to those illustrated in Fig.[Tl which may now be read as follows. The nucleon is lighter than the A 
owing to attraction provided by quark-exchange originating from both scalar- and axial-vector-diquark 
correlations. With increasing quark mass, not only do the constituents of both systems approach 
common masses but the additional attraction is diminished. Hence m^r increases more rapidly than 
does the mass of the Z\, whose simpler structure means it draws attraction from only one source in the 
Faddeev equation. 

There is another interesting consequence of increasing current-quark mass, which concerns the axial- 
vector-diquark content of the nucleon. As Sqq-^-^ becomes smaller, there is less to distinguish between 
the range and nature of the attractive interactions provided by the scalar- and axial-vector-diquark 
correlations. Hence, as shown in Fig.[6j the relative strength of the axial-vector correlation increases 
with m. This relative strength has a material impact on nucleon properties. For example, in the 
electromagnetic form factor calculations described in Ref. [63 | , the photon- nucleon interaction involves 
an axial-vector diquark correlation with 40% probability. This value: is crucial in fixing the location of 
a zero in the ratio Ff''^(Q^/Ff'^(Q^) [H^; determines the x = 1 value for ratio of nucleon structure 
functions F2 /F^ = 0.36 [6^]; and entails that in the nucleon's rest frame just 37% of the total spin of 
the nucleon is contained within components of the Faddeev amplitude which possess zero quark orbital 
angular momentum [29| . 
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Table 4 Row-l: Dressed-quark-core masses for nucleon and A, their first radial excitations (denoted by "*"), 
and the parity-partners of these states, computed with qn = 1.18, qa = 1.56, and the parameter values in 
Eq. (|25p and Table [T] The errors on the masses of the radial excitations indicate the effect of shifting the 
location of the zero according to Eq. pO|l . Row -2: Bare-masses inferred from a coupled-channels analysis at the 
Excited Baryon Analysis Center (EBAC) [63 . EBAC's method does not provide a bare nucleon mass. ROW- 
S'. Bare masses inferred from the coupled-channels analysis described in Ref. [SJ, which describes the Roper 
resonance as dynamically-generated. In both these rows, ". . . " indicates states not found in the analysis. A 
visual comparison of these results is presented in Fig.[71 
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4.4 Radial excitations and parity partners 

In analogy with mesons, the leading Tchebychev moment of the bound-state amplitude for a baryon's 
first radial excitation should possess a single zero. Whilst the truncation we employ cannot generate 
such a zero, it is possible to estimate masses for these states by employing the expedient described 
in Sect. 13.21 With the zero located as prescribed in Eq. (pHl) . no new parameters are introduced. To 
illustrate, the mass of the first radial excitation of the A, m^* , is determined via 

9a -^991+ 1 



^^M^^2^ Jo ('^^("'^^''^991+'"^^*)) ■ (45) 

In a more general setting one might imagine that a baryon's first radial excitation could be an 

admixture of two components: one with a zero in the Faddcev amplitude, describing a radial excitation 

of the quark-diquark system; and the other with a zero in the diquark's Bethe-Salpeter amplitude, 

which represents an internal excitation of the diquark. The procedure in Sect. 13 .2] can conceivably 

distinguish between these components via a mixing term whose strength is oc Eqq_^^ ^qq*^ ■ Owing to 

orthogonality of the two-body ground- and first-radially-excited states, we anticipate that this mixing 

term is negligible. Under this assumption, a baryon's first radial excitation is predominantly a radial 

excitation of the quark-diquark system. Should a state constituted from a radially-excited diquark 

exist, then its mass will be larger because \E^„, /m^„. 1 < \E^„ , /m^„ , 1. 

^ L qq^^ ' qq^^^ ^ qqi+ ' qqi+^ 

Given the preceding discussion, it will not be surprising that we define bound-state equations for 

the parity-partners of the ground- and first-radially-excited-states of the nucleon and Z\-resonance by 

making the replacements 

Eqqjp Eqq^ p , Tllqq^p ^T'qQj^p (46) 

in the appropriate Faddeev equations. For example, we determine the mass of the = | Z\-state 
from 

^ = 17 T'" / da{m 2 + {I - afm^ 3-){am^3- + M)C'^{aA{a,M,rnqq m^s.-)) ■ (47) 

M 'm^q^_ 27r^ Jq ^^i- ^2 ^2 ^2 



4.5 Computed baryon spectrum 

In Table H] we list our computed results for the dressed-quark-core masses of the nucleon and A, 
their first radial excitations (denoted by "*"), and the parity-partners of these states. These masses 
cannot be compared directly with experiment because the kernels employed in their calculation do 
not incorporate the effect of meson loops. However, a fair comparison may be made with bare-masses 
inferred from sophisticated coupled-channels analyses of ttN scattering data up to W < 2 GeV [6^ |63| • 
The predictions of our model for the baryon's dressed-quark-core match the bare- masses determined 
in Ref. |67| with a root-mean-square (rms) relative-error of 10%. Notably, however, we find a quark-core 
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Fig. 7 Visual comparison of our com- 
puted baryon masses (filled circles) with 
bare masses determined in Ref. [63] (filled 
diamonds) and Ref. [63 (filled triangles). 
For the coupled-channels models a sym- 
bol at the lower extremity indicates that 
no associated state is found in the anal- 
ysis, whilst a symbol at the upper ex- 
tremity indicates that the analysis reports 
a dynamically-generated resonance with 
no corresponding bare-baryon state. In 
connection with the i7-baryons the open- 
circles represent our results shifted down- 
wards by 100 MeV. [See discussion after 
Eq. ()48|) .] The fill ed-s quares report masses 
tabulated in Ref. [H. 
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to the Roper resonance, whereas within the Jiilich coupled-channels model this structure in the Pn 
partial wave is unconnected with a bare three-quark state. In connection with EBAC's analysis, our 
predictions for the bare-masses agree within a rms relative-error of 14%. Notably, EBAC does find a 
dressed-quark-core for the Roper resonance, at a mass which agrees with our prediction. 

We also predict dressed-quark-core states associated with radial excitations of the Z\-resonance. 
Allowing for a reduction by < 160 MeV expected from coupled channels effects, our estimate for the 
mass of the three-star zi(1600) Psa-resonance is in agreement with contemporary experiment [42l |. The 
same is true of our result for this state's parity-partner, Zl(1940) D33. 

In this connection we observe that the EBAC analysis does not find any sign of the Zl(1600) P33 
or Z\(1940) -D33 resonances. This is plausibly an indication of a limitation in the method employed to 
complete the difficult task of reaching into the complex plane in order to locate the poles associated 
with these resonances. 

It is also worth remarking that the Jiilich analysis of the / = |-channel has been revisited [68| . 



2 

with new bare masses being reported: 1535 MeV for the Zi(1232) and 3442 MeV for Z\(1700). However, 
this study finds that the bare-mass values depend sensitively upon precisely which channel-couplings 
are included in the model. It is notable that this analysis identifies the Zi(1600) P33 with a broad, 
dynamically-generated resonance. 

With the information now available we can also report dressed-quark-core masses for the decuplet 
f2~ baryon; viz.. 



171(2 ma* 



1.76 2.00 ±0.03 2.56 2.59 ± 0.01 



(48) 



Only four l?-baryons are listed in Ref. f2~(1670), a four-star = |^ state; J7~(2250), a three- 
star state; and i7~(2380), !7~(2470), two two-star states. The spin-parity of the last three resonances 
is currently unknown. To place our computed f2-masscs in context we observe that at m = the 
computed vector-meson dressed-quark-core mass is = 1.13, which is 110 MeV above the experi- 
mental value. Notably, pseudoscalar-meson loop corrections are estimated to reduce the core mass by 
~ 100 MeV [l^;!!^- Furthermore, a similar analysis indicates that, at = 0.5 GeV^, pseudoscalar- 
meson loop corrections in Z\-like systems produce a (—100 MeV) shift in the mass of the baryon's 
drcssed-quark core [ssj . A comparison with experiment is presented in Fig. [7] 



5 Epilogue 

We presented the first DSE-based calculation of the light hadron spectrum that simultaneously cor- 
relates the dressed-quark-core masses of meson and baryon ground- and excited-states within a single 
symmetry-preserving framework. Isospin symmetry was assumed, with m„ = nid = m = 7 MeV pro- 
ducing a physical pion mass; and five parameters were used to define the gap-, Bethe-Salpeter- and 
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Faddeev-equations. In a comparison with relevant quantities, we recorded a value of 13% for the over- 
all root-mean-square- relative-error/degree-of freedom (fms). Notable amongst our results is agreement 
between the computed masses for baryon dressed-quark-cores and the bare masses employed in modern 
dynamical coupled-channels models of pion-nuclcon reactions. 

In connection with mesons we capitalised on recent progress in understanding the far-reaching 
effects of dynamical chiral symmetry breaking within the Bethe-Salpeter kernel to improve upon the 
rainbow-ladder truncation in the scalar and axial-vector channels. This enabled us to obtain TmSqq ~ 
13%, a feature which also has a significant collateral impact on the baryon spectrum owing to the 
connection between the Bethe-Salpeter equations for mesons and diquark-correlations. 

In comparison with relevant quantities our predicted baryon masses yield rmSggg ~ 14%. Fur- 
thermore, our analysis provides insight into numerous aspects of baryon structure. For example, we 
explained that in practical formulations of the baryon bound-state problem, there are two distinct 
types of pseudoscalar-meson-loop correction: one intrinsic to the gap equation; and another restricted 
to bound-state kernels. It is the latter which is expressed in the formulae typically used to estimate 
meson-cloud contributions to baryon masses. 

We also demonstrated relationships between the masses of the nucleon and A (m^r, m^), and 
those of the dressed-quark and diquark correlations they contain (M, rnqq^^ , m^q^^ ). For example, we 
established a causal connection between [m^ —misi] and [mqq_^^ ~'^99o+]' demonstrated that ttln ^ 3Af 
and TO/i AI + 'mqq_^^ ; illustrated the simplicity of the Z\'s internal structure and the consequences of 
this; and showed that the relative strength of the axial-vector diquark-correlation within the nucleon 
grows with current-quark mass. 

At the core of our analysis is a symmetry-preserving treatment of a vector- vector contact interaction. 
Our body of results confirms that this is a useful tool for the study of phenomena characterised by 
probe momenta less-than the dressed-quark mass. It is now important to use this foundation in the 
computation of baryon elastic and transition form factors. That will provide information which is crucial 
in using experimental data on such observables as a tool for charting the nature of the quark-quark 
interaction at long-range [T^. 
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A Euclidean Conventions 

In our Euclidean formulation: 

4 

p-q = '^Piqi; (A.i) 

1=1 

{Tm. 7^} = 2 (5,1^ ; 7^ = 7,1 ; ct,,^ = ^ [7^, 7^] ; tr [757M7f 7p7<t] = -4 e^-^pa , £1234 = 1 . (A.2) 
A positive energy spinor satisfies 

u{P, s) [i-y- P + M) = ^ {i-y- P + M) u{P, s) , (A. 3) 

where s = ± is the spin label. It is normalised: 

m(P, s) u{P, s) = 2M , (A.4) 

and may be expressed explicitly: 



u{P, s) = VAf - i£ o- P , (A.5) 



with £ = iVP^+AP, 

X+ = (;).X- = (?). (A.6) 
For the free-particle spinor, s) = u{P, 3)^74. 
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The spinor can be used to construct a positive energy projection operator: 

A negative energy spinor satisfies 

v{P, s) [i-y-P- M) = = (i7 ■ P - M) v{P, s) , (A.8) 

and possesses properties and satisfies constraints obtained via obvious analogy witli u{P,s). 
A cliarge-conjugated Bethe-Salpeter amplitude is obtained via 

r{k;P)^C'^r{-k-PfC, (A.9) 

where "T" denotes a transposing of all matrix indices and C = 7274 is the charge conjugation matrix, — —C. 
We note that 

CSJC = -7m, [C,75l =0. (A.IO) 

In describing the A resonance we employ a Rarita-Schwinger spinor to unambiguously represent a covariant 
spin-3/2 field. The positive energy spinor is defined by the following equations: 

(i7-P + M)u^(P;r) = 0, -f^u^{P;r) = , P^,u^{P-r) = Q, (A.ll) 

where r = —3/2, —1/2, 1/2, 3/2. It is normalised: 

u^,{P-r')u^{P-r) = 2M , (A.12) 

and satisfies a completeness relation 

3/2 

2M 



0/ z 

^ ^ u^{P-r)u,{P-r) = A+{P)R^,, (A.13) 

r = -3/2 

where 

1 2 - - 1 - 

Pm" = S^^I-D - g7M7^ + -^Pi^PulYy - ig[PM7^ " -P.' 7m] . (A.14) 

with P'^ = —1, which is very useful in simplifying the positive energy Z\'s Faddeev equation. 



B Bethe-Salpeter equations 

B.l Pseudoscalar mesons and scalar diquarks 
The explicit form of Eq. (fT2)) is: 







1 




i^7r 

r^EF 
















r^EF . 






(■P). 



where 



^EE 

K-EF 

K-FE 
fCpF 



da C"'{uj{M'^,a, P^)) - 2q(1 - a) P^ C^{lu{M^, a, P^)) 



= P' 
~ — 2/Cfb 



daCi {uj{lVr,a,P'^)), 



daCTiuj{M'^,a,P^)), 



(B.l) 

(B.2) 

(B.3) 

(B.4) 
(B.5) 



with Ci{z) =Ci{z)/z. We used Eq. (|22l) to arrive at this form of ICff- It follows immediately that the explicit 
form of Eq. (fT4)l is: 



P 



I^EE i^EF 
'^FE I^FF 



E,,„+ (P) 



P 



(B.6) 



Equations ([B.ip and (|B.6|) are eigenvalue problems: they each have a solution at a single value of P^ < 0, at 
which point the eigenvector describes the on-shell Bethe-Salpeter amplitude. 
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In the computation of observables, one must use the canonically-normalised Bethe-Salpeter amplitude; i.e., 
Ftt is rescaled so that 

Pa = N,tr f^r^[-P)^S{q + P)r^{P)S{q), (B.7) 



where Nc = 3. For the pion in the chiral limit, this expression assumes a particularly simple form; viz. 



(B.l 



The canonical normahsation condition for the scalar diquark is almost identical to Eq. (|B.7|) . with the only 
difference being the replacement Nc = 3 — >■ 2. 

In order to estimate the mass of the first radial excitation of the pion, we use the following kernel 



f^EE 
K-EF 
iCpE 



da 



P''{uj{M^,a, P^)) - 2a(l - a) P^ J^^{uj{M^, a, P')) 



(a;(Mo^a,P')), 



(B.9) 

(B.IO) 

(B.U) 
(B.12) 



It is conceived so that the leptonic decay constant of the radially-excited pseudoscalar meson vanishes in 
the chiral limit, which is a consequence of the axial- vector Ward-Takashi i dent ity [5^ . The kernel for the 
radially-excited scalar diquark is obtained through obvious analogy with Eq. (|B.6|I . 



B.2 Mesons and diquarks with J = 1 

In the treatment of Eq. Q using the rainbow-ladder truncation the vector-meson Bethe-Salpeter amplitude 
has a particulary simple form; viz., 

rr (P) = 7^ -El- (P) , = 0. (B.13) 

Hence the explicit form of Eq. ([T]) for the ground-state vector-meson, whose solution yields its mass-squared, is 

l + K''{-ml-)=0, K''{P^) = -^—^ [ daa{l-a)P^cT{io{M\a,P^)). (B.14) 

STT^m^G Jo 

Equation (|22p was used to express the BSE in this form. The BSE for the axial-vector diquark again follows 
immediately; viz., 

l + i^"(-"4,+ )=0. (B.15) 
The canonical normalisation conditions are readily expressed; viz.. 



_1 



2 ="9mS7B2^'(^ 



772 ^ dP'^ 



(B.16) 



We emphasise that the vector-meson and axial-vector-diquark BSEs only assume such particularly simple 
forms in the rainbow-ladder truncation. Even with a momentum-independent interaction, vector meson and 
axial-vector diquark Bethe-Salpeter amplitudes possess two Dirac covariants immediately upon inclusion of 
next-to-leading-order corrections to the quark-gluon vertex; viz., 

r!.' (P) = 7m iSi- (P) 7m -El- (P) + ij^a^.P.F,- (P) , 7^P^ = 0. (B.17) 

Similar observations hold for a g^D{p — q) ~ (5*(p — q) interaction 0; [Toi: [Tll|. 

Again owing to the simplicity of the interaction, the Bethe-Salpeter amplitude for an axial- vector meson is 

/^m*(^')= 757m-Bi+(P). (B.18) 
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In this case dressing the vertex does not generate new covariants because a momentum-independent intera ction 
cannot generate a Bethe-Salpeter amphtude that depends on the relative momentum. Inserting Eq. HB.18|I into 
Eq. @ yields the following BSE: 



1 + A-"! (-m?+ ) = , (P^) = 
It follows that the vector-diquark mass is determined by 

The canonical normalisation conditions are 

1 o 2 d 

— — 9rr 

^1+ 



P2 = - 



(B.19) 



(B.20) 



(B.21) 



B.3 Scalar mesons and pscudoscalar diquarks 

The Bethe-Salpeter amplitude for a scalar meson is 

roAP) = lBEo+(P). (B.22) 

As with axial-vector mesons, dressing the vertex does not generate new covariants. Inserting Eq. (|B.22[I into 
Eq. ((31 yields the following BSE: 

16 1 /"A q^ + q.P-M^ ^^24) 



3 m% J (27r)4 [{q + P)2 + AP][q^ + M^] ' 

Now consider Eq. p8|l : if one sets P ^ = — 4M^ in that chiral limit identity, then one finds after just two lines of 
algebra that it is equ ivalent to Eq. (|B.24[I . Hence, for m = the treatment of Eq. ((TJ using the rainbow- ladder 
truncation yields [33| 

mo+ = 2 M . (B.25) 

For ge neral values of the current-quark mass, using our symmetry-preserving regularisation prescription, 
Eq. (|B.23I) can be written 

l + K^-ml+) = 0, K''{P^) = -^ [ da [c'"MM^ «, P^)) - 2C^MM^ a, P^))] . (B.26) 



'G Jo 

It follows that in the rainbow-ladder truncation the mass of a pseudoscalar diquark is determined by 



l + iA-(_^2^ ) = 0. (B.27) 



The canonical normalisation conditions are 

1 9 2 d ,^..„2 



El, 2'^^dp-^^'^^^ 



0+ " 



(B.28) 

P2 = -mL 



C Faddeev Equation 

C.l General structure 

The nucleon is represented by a Faddeev amplitude 

^ = <I'i + <I'2 + <I'3 , (C.l) 

where the subscript identifies the bystander quark and, e.g., !Z'i,2 are obtained from by a cyclic permutation 
of all the quark labels. We employ the simplest realistic representation of •f'. The spin- and isospin-1/2 nucleon 
is a sum of scalar and axial-vector diquark correlations: 



(C.2) 
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with {pi,ai,Ti) the momentum, spin and isospin labels of the quarks constituting the bound state, and P = 
Pi + P2 + P3 the system's total momentum. 

It is conceivable that pseudoscalar and vector diquarks could play a role in the ground-state nucleon's 
Faddeev amplitude. However, they have parity opposite to that of the nucleon and hence can only appear in 
concert with nonzero quark angular momentum. Since one expects the ground-state nucleon to possess the 
minimum possible amount of quark orbital angular momentum and these diquark correlations are significantly 
more massive than the scalar and axial- vector (Table [S]), they can safely be ignored in computing properties 
of the ground state. 

The scalar diquark piece in Eq. (|C.2|) is 



AAf (p„a„rO = [r°^(ip[i2];if)]?;4 A°\k) [S{e-P)uiP)]Z , (C.3) 

where: the spinor satisfies (App.|X|) 

{i'y- P + M) u{P) = = u{P) [i'y- P + M), (C.4) 

with M the mass obtained by solving the Faddeev equation, and it is also a spinor in isospin space with 
Lp+ = col(l,0) for the proton and — col(0, 1) for the neutron; K — pi + p2 =■ P{i2}, P[i2] = Pi — P2, 
i ■■= (-P{i2} + 2p3)/3; 

is a propagator for the scalar diquark formed from quarks 1 and 2, with mg+ the mass-scale associated with 

this correlation, and P" i s the canonically-normalised Bethe-Salpeter amplitude describing their relative mo- 
mentum correlation. Sect. IB. l1 and 5, a 4 x 4 Dirac matrix, describes the relative quark-diquark momentum 

correlation. The colour antisymmetry of ^f'a is implicit in P'^ , with the Levi-Civita tensor, £0^0203, expressed 
via the antisymmetric Gell-Mann matrices; viz., defining 

{H^ = -iA^i^^ = ^A2}, thenecic^cs = iH''')c,c2- (C.6) 

The axial- vector component in Eq. (|C.2|I is 



M'\p,,a,,n) = [t'r^''{^Pii2y,K)]l\X ^tiK) [^.{t, P)u{P)]l% , (C.7) 



where the symmetric isospin-triplet matrices are 

t = ^{-r +T ), t =T , t = — (r - r ) , (C.8) 
and the other elements in Eq. (|C.7|) are straightforward generalisations of those in Eq. (|C.3|) with, e.g., 

< (K) = ^2^^ 2 (s.-^ + ^] ■ (C.9) 

Since it is not possible to combine an isospin-0 diquark with an isospin-1/2 quark to obtain isospin-3/2, 
the spin- and isospin-3/2 A contains only an axial- vector diquark component 

>l>f{p,,a,,n)=Vl+. (C.IO) 

Understanding the structure of the A is plainly far simpler than in the case of the nucleon since, whilst the 
general form of the Faddeev amplitude for a spin- and isospin-3/2 can be complicated, isospin symmetry means 
that one can focus on the A'^'^ , with its simple flavour structure, because all the charge states are degenerate: 

= [t+r',\lpii2y,K)]l\X AIUk) [V,,{e;P)u,{P)^+]ll , (C.ll) 

where Up[P) is a Rarita-Schwinger spinor, Eq. ()A.11|) . 

The general forms of the matrices S(£;P), Al^ii'jP) and T>vp{l\ P), which describe the momentum-s pace 
correlation between the quark and diquark in the nucleon and A, respectively, are described in Refs. [29l: I23||. 
The requirement that S{1; P) represent a positive energy nucleon entails 

S{l;P)=si{l-P)I^ + {i'yi-i-PIj,) S2{t,P), (C.12) 

where (/D)rs ~ Srs, = 1, P^ — —1. In the nucleon rest frame, Si,2 describe, respectively, the upper, lower 
component of the bound-state nucleon's spinor. Placing the same constraint on the axial-vector component, 
one has 

6 

n^Y. P-(^' ^5 P) , i = +, 0, - , (C.13) 
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where 



= 7 • P,, Al = ~iP, , Al=-fi^ 
At =iii, Al^^^-Al,Al= i7^7 ■ 



■At. 



Finally, requiring also that 'D^p{£;P) be an eigenfunction of yl+(P), one obtains 

O.p {£; P) = 5"" {£; P) <5.p + j^A^ {£; P) if , 

with 5^ and A^ given by obvious analogues of Eqs. (|C.12|) and (|C.13[) . respectively. 
One can now write the Faddeev equation satisfied by 'I'a as 



S{k-P)u{P) 
Al(k-P)u{P) 



= -4 



M{k,£;P) 



S{£;P) u{P) 
Aii£;P)uiP) 



The kernel in Eq. IfOJE} is 



Mik,£;P) 



Moo (Moi) 



with 



Moo = r°*{k, - £,J2; £,,) S'^ii,, - k,) r°*{£g - k,j2; -fc,,) S{£,) {£,,) , 



where: £g = £ + P/3, kg = fc + P/3, £gg = -^ + 2P/3, kgg = -fc + 2P/3 and the superscript "T" denotes 
transpose; and 

{M0l)i = rl^{kg - £gg/2- £gg)S^{£gg " kg) r°^{£g - kgg / 2' -kgg) S {£ g) A^^, {£ gg) , 

{Ml0)i = r°^{kg - £gg/2; £gg) ^ kg) F^'^ {£ g - kgg/2; ~ kgg) S{£g) Zi^^ (£gg) , 

{Mll)j,y — "t'' Pp (kg — £gg/2; £gg) S {£gg — kg) t (£g — kgg/2; ^kgg) S(£g) A (£gg) ■ 

The A's Faddeev equation is 

d^£ 



Vxpik;P)up{P) = 4 



(27r) 



jMt^k,£;P)V^^{£;P)u^{P)., 



with 



A4a/j, — "t J^a {^q ^qq /^'i ^qq) S {^qq kg) t /\ [£q kqq/2; kqq) S {£q) ^J.{^qq) ■ 



(C.14) 
(C.15) 

(C.16) 

(C.17) 

(C.18) 
matrix 

(C.19) 
(C.20) 
(C.21) 

(C.22) 
(C.23) 



C.2 Explicit form of the nucleon's Faddeev equation 

Using Eq. (|31|l . the nucleon's Faddeev amplitude simplifies and can be written in terms of, Eqs. (|C.12|I . 

5(P) = s{P) 7d , ^;(P) = al(P)757M + 4(P)75Pm ,^ = +,0 . 
The mass of the ground-state nucleon is then determined by a 5 x 5 matrix Faddeev equation; viz., 



/ siP) \ 


( 






1i 




^01 


\ 


/ 4P) \ 


a+(P) 




-V2Al«, 





\/2 J^aiai 





^2 A^^j,2 




a+iP) 


aliP) -- 




^10 




A'" 

aiai 




A" 




aliP) 


4iP) 














\/2 A^2a2 




4iP) 


V aUP) J 


V 




\/2 Kllai 




\/2 Aa2a2 


A" 


/ 


V aUP) J 



where: cjv ~ g^/(47r^Af), 

o-jv = (^N{a, M, mgg^^ ,mjv) := (1 — a) + am^g^^ 
and 

A° 

Abb = cjvP^ 



i{l - a)m% , (t]v = ON{a, M, 



,niN) 



^00 , ^00 , T^OO 



EF 
1 



11o+ 



daCi{aj^){amN + M) . 



Hem . 

(C.24) 



(C.25) 



(C.26) 

(C.27) 
(C.28) 
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K 



(C.29) 



K°A = cnF^,^^ y daCiia%){amM + M) ; (C.30) 

-^sai — ^SEai + -^sj;-ai : (C.31) 



Enn..^ En 



Kf^a, = cjv^^ai-^^ / daCi(a]v)(m^, , (3M + amiv) + 2a(l - afm%) , (C.32) 

"^991+ Jo 
01 ^<i<io+^<iQi+ rriN 



^ -Civ % irf / (M + 3amjv) + 2(l-a)"A/m^); (C.33) 

^si, = ^"^a, + ^s^a, , (C.34) 
»1 



°^ - cjv— ^^2± — / daCi((T]v)(am]v - Af)((l - a)^m?^ - m^,^^ ) , (C.35) 



'991+ Jo 



J^'^paa = cjv-^^2V"^^"T7" / c'"Ci(o-jv)(l - Q!)(amjv - Af)((l - a)^mir - m^„ ; (C.36) 



-^ais — -^aiSE +-^aiSF j (C.37) 



^10^^ = ^:E!!2±f!£l± ^ daCi(a^)(am,v + M)(2m^,^^ +(l-afm^), (C.38) 



'991+ "'0 

_ CJV -^99o+ -^99i+ mjv 
is^ - ^^1^-^ 



dQCi(o-SV)(l - a){2mlq^_^ + (1 - a)^mif)(am]v + M) ; (C.39) 



Kll = i^I^.E + , (C.40) 

»i 



= ^ °, / dQCi(a^)(amiv + Af)K,^^ _4(l-a)^m^), (C.41) 



KlU = ^^2!2±f!£l±!^ f daC,{a%){l-a){^ml,^^ - 2{l - afml){am^ + M) ; (C.42) 

".1=-^^^/ daCi(a^)[2m^,^^(Af-amiv) + (l-a)2m^(M + 5amjv)l; (C.43) 

•J '>TI,qq-^+ Jo 
11 ZCJV -^^991+ / j_ — 



a^a, = / daCi(a^)(-m^,^^ + (1 - a^)m^)(amiv - M) ; (C.44) 

J ^^991+ Jo 



11 _ CN -^99i + 



a2ai o 2 

O lll-qq-^^ Jq 



daCi(al,)[mlq^^ (llamjv + M) - 2(1 - afmi^{7amN + 2M)] ; (C.45) 



5CJV -^99i + 



3 "^99?^ ^0 



daCi{aff){mqq^^ — (1 — a) m,^){amN — M) . (C.46) 

^+ 

This kernel was computed following t he p rocedure detailed for the Z\-resonance in Sect . 14.11 During this process 
we employed the replacements in Eq. (|40|l . their analogues involving the scalar-diquark's momentum, A'o+, and 
KqJt ■ K^Jr — >■ (1 — af'P^ . In the present context, of course, = —m%. 

Given the structure of the kernel, it is not surprising that the eigenvectors exhibit the pattern 

a+ = -^2a°, i = 1,2. (C.47) 

For example, at the mass presented in Table 21 the nucleon's unit-normalised Faddeev amplitude is 

if a? a2 
0.88 0.38 -0.27 -0.065 0.046 



S flj fli flj "2 ^(J ^g-J 



The axial-vector-diquark correlation provides 22% of the unit normalisation. This is discussed further in con- 
nection with Fig.[6l 
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